Apical oxygens and correlation strength in electron and hole doped copper oxides 
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We use the Local Density Approximation in combination with the Dynamical Mean Field Theory 
to carry out a comparative investigation of a typical electron doped and a typical hole doped 
copper oxide, NCCO and LSCO respectively. The parent compounds of both materials are strongly 
correlated electron systems in the vicinity of the metal to charge transfer insulator transition. In 
NCCO the magnetic long range order is essential to open a charge transfer gap, while Mott physics 
is responsible for the gap in LSCO. We highlights the role of the apical oxygens in determining 
the strength of the correlations and obtaining overall good agreement between theory and several 
experimentally determined quantities. Results for optical conductivity, polarized X-ray absorption 
and angle resolved photoemission are presented and compared with experiments. 



INTRODUCTION 

Since their discovery, the electronic structure of the 
high temperature superconductors has been a subject of 
intensive theoretical attention as well as controversy, a 
situation that continues even today. A landmark ques- 
tion to understand these materials is how their physical 
properties follow from their electronic structure and to 
which extent simplified descriptions in the form of model 
Hamiltonians describe their basic physical properties in 
the normal state. 

It is generally accepted that the physics of the copper 
oxide based high temperature superconductor families is 
captured by the copper-oxygen layers and the relevant 
degrees of freedom are the copper dx2_y2 orbitals and 
the oxygen p^; and Py orbitals [T, ^ . Numerous studies 
have demonstrated that this model captures qualitatively 
the physics of the copper oxide planes. However, there 
have been several proposals that the p±z oxygen orbitals 
and the dsz2-r2 copper orbital play an important role in 
the onset of orbital current order [3 or for the existence 
of superconductivity [4, 5 . 

Seeking a simpler low energy description, several stud- 
ies [6l [7] have shown that the Hubbard model describes 
some qualitative properties of the copper oxygen layers. 
However, the precise energy range over which the de- 
scription is valid and the quality of this description for 
different physical observables, is still a subject of active 
research. 

A related question is hence how to map the copper 
oxide layers onto the various effective Hamiltonians and 
what is the effective strength of the Coulomb correlations 
in these systems. There are various approaches tack- 
ling this issue, ranging from ab-initio methods to model 
Hamiltonian studies. 

The theoretical studies usually consist of a two stage 
process: In the first step ab-initio approaches use an ap- 
proximate technique, such as constrained DFT in the 
LDA or GGA approximation, or quantum chemical meth- 
ods, to derive the parameters and the form of the effective 
Hamiltonian JSHTO]. In the second step, model Hamilto- 



nian based approaches compute various observables in 
the framework of a given model Hamiltonian and given 
approximation technique. The parameters entering the 
model Hamiltonians are determined by comparing the 
results side by side with experiment. 

Numerous research efforts notwithstanding, even the 
basic question of the strength of the correlations in the 
copper oxide materials and the origin of the insulating 
gap in their parent compound has not been fully eluci- 
dated. There are two opposing physical pictures describ- 
ing the origin of the insulating gap in these materials. 
In the so called Slater picture, the insulating behavior 
is understood as the result of a doubling of the unit cell 
caused by antiferromagnetic long range order. In the so 
called Mott picture, the insulating behavior is the result 
of the local blocking of the electron propagation due to 
the strong Coulomb repulsion. In the latter picture, the 
insulating behavior originates from the localization of the 
electron and is not tied to any specific form of magnetic 
long range order. Hence antiferromangetic long range or- 
der arises as a secondary instability. In the presence of 
magnetic long range order, the unit cell is doubled and 
the two pictures. Slater and Mott, are continuously con- 
nected: No physical observable can provide a sharp dis- 
tinction between the two, both magnetic order and block- 
ing contribute to the insulating behavior. Mean field the- 
ory treatments allow to study the paramagnetic state as 
an underlying normal state mean field solution, which 
can supports a sharp transition between a paramagnetic 
metal and a paramagnetic insulator. This solution is not 
realized when other more stable mean field solutions sup- 
porting long range order exist, but within a mean field 
framework, it can still be used to draw a sharp distinc- 
tion between Slater insulators and Mott insulators, by 
investigating if the ordered state is derived from a metal- 
lic or an insulting paramagnetic solution. The pioneer- 
ing work of Zaanen Sawatzky Allen [11] and their sharp 
(so called ZSA) boundary between charge transfer metals 
and charge transfer insulators can be viewed in this light. 

A simple argument can be formulated for the Hub- 
bard model to estimate the strength of the correlation 
of cuprates. Within a one band Hubbard description. 



parametrized by a bandwidth W ^ 8t and a Coulomb 
repulsion U, the insulating gap of the paramagnetic in- 
sulator is U — W and the super-exchange is J = At^ /U. 
For cuprates, the gap is around ~ 2eV and J ~ O.leV, 
and therefore it is found that the repulsion is of the order 
U /W ~ 1.5, which is above but not far from the Mott 
boundary U ^W^ and hence the cuprates are in a regime 
of intermediate correlation strength. Conclusions on the 
placement of cuprates in a regime of intermediate cor- 
relation strength were also reached by numerical studies 
(for a review see for example Ref. [12j. 

The strength of the correlations was also studied in 
the three band theory. In particular, large N slave bo- 
son mean field theory of a three band model (with no 
oxygen-oxygen transfer integrals) [13] of the copper ox- 
ides obtained a sharp transition between the metal and 
the charge transfer insulator in the paramagnetic phase. 
This metal-to-charge-transfer-insulator transition paral- 
lels the Brinkman Rice transition in the Hubbard model. 
It was found that there is indeed a correspondence be- 
tween the two critical Coulomb repulsions of the Hubbard 
model Uci (the minimal Coulomb repulsion that supports 
the paramagnetic insulator) and Uc2 (the minimal repul- 
sion that does not support a paramagnetic metal) with 
the two critical charge transfer energies Ad and Ac2, 
where the paramagnetic insulating state and the para- 
magnetic metallic state are destroyed, respectively. The 
critical value of the charge transfer energy in the three 
band theory (Ac2) plays the role of the critical U {Uc2) 
of the Hubbard model [13] . Although a strong particle- 
hole asymmetry is expected in the three band theories, 
since doped electrons reside mainly on copper sites while 
doped holes reside mainly on the oxygen sites, it was 
shown that the resulting quasiparticle band structure was 
surprisingly particle- hole symmetric. Indeed, it is due to 
the strong copper-oxygen hybridization which results in 
the formation of Zhang-Rice singlets and to the quasi- 
particles that involves copper spin and oxygen charge. 

Other more realistic treatments of the three band de- 
scription of the copper oxides were carried out within the 
slave boson framework. Some include the oxygen dis- 
persion [14] [15], additional copper and oxygen orbitals 
[T6] , short range magnetic correlations [TTl [18] , the near- 
est neighbor Coulomb interactions [19] or the electron 
phonon coupling [20 . Within slave boson mean field the- 
ory of the three band model, the parent compound of hole 
doped cuprates LSCO was located close but above (i.e. 
on the insulating side of) the metal to charge transfer 
insulator transition boundary [21 . 

The development of Dynamical Mean Field Theory 
[22] and their extensions [23H25] opened new avenues to 
advance our qualitative understanding of the electronic 
structure of the cuprates and its quantitative description. 
DMFT goes beyond slave bosons theories: This method 
treats both coherent and incoherent features on the same 
footing, whereas slave bosons theories are not able to 



capture the coherent character of the physical solution. 
DMFT successfully describes the Mott transition of the 
Hubbard model and gives a deeper understanding of the 
Brinkman Rice transition [22j [26]. Single site DMFT 
calculations can also be extended to more precise cluster 
calculations. In particular, cluster corrections in DMFT 
allow to assess the validity of the single site calculations. 

The question of the strength of the correlations was 
also addressed by DMFT studies of multi-band model 
Hamiltonians for the copper oxides planes [271 [2Zl 121] • 
The phase diagram, with respect to the charge transfer 
energy and the Coulomb repulsion of the copper orbitals 
was studied. The boundary between the metallic and the 
charge-transfer-insulator solutions was located, as well as 
the crossover line between the charge-transfer-insulator 
and the Mott insulator. A full phase diagram of a cop- 
per oxide model, and a very complete analogy with the 
DMFT studies of the Mott transition in the Hubbard 
model was also performed recently [29] , 

Combination of DMFT with electronic structure meth- 
ods, such as LDA, allow to combine the ab-initio and 
the model Hamiltonian viewpoint in the LDA+DMFT 
framework. The LDA+DMFT method [23| allow in par- 
ticular to determine the strength of correlations for spe- 
cific materials, like NCCO and LSCO, and there are still 
several important issues unresolved regarding how this 
picture is connected to the cuprates: a) How should the 
different materials be placed in the qualitative ZSA phase 
diagram. Should the parent compounds of the copper ox- 
ide materials be thought as Slater or Mott /charge trans- 
fer Insulators, b) What significant differences are there 
in the different level of description, mainly what are the 
differences between the one band and three band the- 
ory, c) What is the quality of the description of the vari- 
ous experimental observables, for the different low energy 
models involving a different number of bands, and finally 
d) Can one obtain a consistent picture of the spectro- 
scopies of hole and electron doped cuprates using a first 
principles method. 

In regards to point a), the issue is still under debate. 
Using the analysis of model Hamiltonians, Refs ^29i |30] 
classify the parent compounds of electron and hole doped 
compounds as Slater insulators in the metallic side of 
the ZSA phase diagrams. Previous LDA+DMFT studies 
classified LSCO [31] as a Mott insulator (or more pre- 
cisely as a charge transfer insulator), and NCCO [32 was 
identified as Slater insulator. On the other hand, the first 
principles study of [331 found NCCO's parent compound 
to be Mott insulator. Finally, the phenomenological anal- 
ysis of experimental data in Ref [34 concludes that the 
analysis of the optical data requires comparable strength 
of interactions in hole and electron doped cuprates. 

In regards to point b) the energy range over which 
the spectral functions of the three band model are repro- 
duced by the one band model is also subject of contro- 
versy. For example, Refs. [29l [30] assert that the three 



band model and a one band model of the copper oxides 
are quantitatively equivalent up to a scale of 4eV, while a 
similar analysis by A. Macridin et al. concludes that the 
range of validity of the Hubbard model is much smaller, 
and is of the order of 0.5 eV [35 . An additional contro- 
versy, regarding the differences between the one band and 
the three band theory, is connected to the values of the 
insulting gap in the paramagnetic insulator. In particu- 
lar, in a one band theory near Uc2^ the gap of the para- 
magnetic solution is substantially smaller then the gap 
of the antiferromagnetic solution for the same parame- 
ters [30[ [Ml El] . This is not the case in the three band 
model where antiferromagnetism increases the value of 
the paramagnetic insulating gap by less than 15% [311. 
The latter statement is controversial with respect to Ref 
p9^ , which argued that the one band model and the three 
band model are similar in their physics up to energies as 
large as twice the gap. 

Regarding point c), many authors used the theoretical 
results obtained with the Hubbard model to fit exper- 
iments. For instance Ref. [38] considered a treatment 
of the Hubbard model with variational cluster pertur- 
bation theory (VCPT). It could match experiments but 
it required a sensible dependence of the Hubbard U on 
the level of doping. On the other hand, an approximate 
diagrammatic treatment of the one band model [39 indi- 
cates that the experimental optical spectra and the dis- 
persion gni [41] can be reproduced without a doping de- 
pendent U. 

Note that a three band LDA+DMFT study was able 
to describe the experimental data of NCCO without hav- 
ing to invoke doping dependent parameters [32 . More- 
over, using the same technique, a successful description 
of both the integrated optical weights and the magnitude 
of the optical conductivity below the charge transfer gap 
for LSCO and NCCO was obtained [31]. This approach 
however does not give the correct magnitude of the opti- 
cal conductivity of LSCO for energies of the order of the 
charge transfer gap, suggesting that additional orbit als 
might play an important role in LSCO. 

In this paper we reconsider these issues building on 
our earlier work of Refs [31] and [32]. We use an ab- 
initio approach, e.g. the Local Density Approxima- 
tion combined with the Dynamical Mean Field Theory 
(LDA+DMFT) [23] to study the electronic structure of 
NCCO and LSCO. The comparative study of two typical 
cuprate compounds LSCO and NCCO allow us to place 
firm bounds ascertaining the importance of correlations 
in the cuprates. The good agreement between theory 
and experiments, achieved using single site LDA+DMFT 
within a multiband framework, for two different com- 
pounds, is a significant results and illustrates the power 
of this new first principles approach to correlated materi- 
als. Our results contrast with a recent phenomenological 
analysis of optical data, of electron and doped cuprates: 
They concluded that for a one band theory, vertex correc- 



tions beyond single site DMFT were required to obtain 
a reasonable fit to experiments [34:\ . 

In this work, we include also the apical oxygens and 
an additional copper orbital {dsz2-r2 ), which were not 
included in our previous study of LSCO [31]. This ex- 
tends the quantitative agreement between theory and 
experiments to a broader energy range. It sheds light 
on why the three band model description of NCCO [32j 
agrees with experiments up to larger energy scales than 
for LSCO [31]. NCCO, in the T' structure, lacks api- 
cal oxygens and is therefore well described by the three 
band model up to a much higher energy scale, justifying 
the excellent agreement between theory and experiment 
found in Ref [32j . 

We confirm that for integration cutoffs smaller than 
half the gap of the parent compound, the additional api- 
cal oxygen degrees of freedom do not affect the integrated 
optical weights in LSCO. This validates the analysis car- 
ried out in Ref [[31] to extract the strength of correlations 
in these materials, which was based on analysis up to an 
energy scale of half the charge transfer gap, i.e. 1 eV. The 
apical oxygens however substantially affects the shape of 
the optical conductivity and it strongly modifies the op- 
tical conductivity of LSCO around 2 eV. While the hole 
occupancy of the Pz orbital is small, it has a clear effect on 
the optical conductivity of the hole doped cuprates. We 
show here that within the first principles LDA+DMFT 
framework, extensions of the model to include further 
orbitals or longer range correlations (using cluster exten- 
sion) consistently improves precision of the calculation 
and improves agreement with experiments. 

The organization of our paper is the following. In sec- 
tion 1 we describe the phase diagram of NCCO and LSCO 
within single site and two site cluster DMFT, and high- 
light the role of magnetic order and singlet correlations 
in these materials. In section 2 we present the evolution 
of the angle resolved photoemission spectra (ARPES) of 
these materials, stressing the various features of the the- 
ory that require a description beyond static mean field 
theory and its comparison with the experiments. Details 
regarding the formula used in section 2 to compare the 
spectral weight with experiments are given in Appendix 
A. Section 3 focuses on the optics and how the evolution 
of the optical properties with doping and temperature 
in NCCO and LSCO can be understood as the result of 
placing the two compound on two different sides of the 
ZSA boundary, once the band structure of both materials 
is taken into account. The inclusion of the apical oxy- 
gens in the theoretical modeling is important to obtain 
correct results for various physical quantities, including 
the shape of the optical conductivity and the integrated 
optical spectral weights with a cutoff of the order of the 
charge transfer gap. We also discuss the connection of the 
optical conductivity with various features in the ARPES 
spectra, and compute the occupation of the different or- 
bitals which are relevant to the XAS spectra. There is 



one parameter, the double counting correction, whose de- 
termination in the LDA+DMFT approach is not unique. 
We thus examined in the last section the dependence of 
our results on this parameter and used this dependence 
to estimate the proximity of both NCCO and LSCO to 
the ZSA boundary. The tight binding parametrization 
obtained makes contact with model Hamiltonian studies 
and is reported in Appendix B. We conclude with some 
outlook for further work. 



FORMALISM 

LDA+DMFT uses first principles density functional 
theory methods to extract the hopping parameters of 
the model, which is subsequently solved using DMFT 
and its extensions. The LDA calculation was carried 
out with the PWSCF package [42 , which employs a 
plane- wave basis set and ultrasoft pseudopotentials [43] . 
Downfolding to a three band model, containing cop- 
per dx2-y2 and two oxygen p^r orbit als was performed 
by the maximally localized Wannier functions (MLWF) 
method [44, ,45^ . The downfolded LDA band structure of 
Nd2Cu04 (NCCO) and La2Cu04 (LSCO) (see Appendix 
B) results in the following three band Hamiltonian: 



nt = 



tTci_c,^, 



iaae{pxjPy) 



{ed-E^')Y,n^da (1) 



where i and j label the Cu02 unit cells of the lattice, 
and t^j are the hopping matrix elements, e^ and Cp are 
the on-site energies of the d and p orbitals, respectively. 
Finally, we note that the charge transfer energy e^ — Cp 
plays the same role as U in the single band Hubbard 
model, as seen for example in slave bosons mean-field 
studies [2T. 

To this Hamiltonian, we add the onsite Coulomb re- 
pulsion U on the dx2-y2 orbital 



^u = ^rf ^ riid^hidi 



(2) 



where the value of Ud = 8eV. The LDA+DMFT method, 
accounts for the correlations which are included in both 
LDA and DMFT by a double counting correction to the 
(i-orbital energy, Edc = Udijid — 0.5), which amounts to 
a shift of the relative positions of the d and p orbitals. 
Here we take Ud to be the occupancy of the correlated 
orbital in the parent compound, which gives the double 
counting corrections Edc = 4.8eV (3.12eV) for NCCO 
(LSCO). 

The LDA downfolded parameters are shown in table [T| 
which we find to be close to those extracted by other first 
principles methods. An extended 6-band model, that 



TABLE L LDA band calculations gives us different set of 
parameters for LSCO ,NCCO and PCCO compounds. The f 
states of the Nd and Pr atoms have been treated as core states 
and are not treated as valence state. Long-range hopping (not 
shown) are also considered within the calculations. The am- 
plitude of the nearest neighbors hoppings (tpp^tdp), the LDA 
on-site energies (cp^e^) and the on-site repulsion Ud are shown 
in this table in eV. In electron notations the bonding orbitals 
enter the Hamiltonian with a negative transfer integral sign, 
and the anti-bonding orbitals with a positive sign. 



Compound 


e'd - 4 [eV] 


Ud [eV] 


td, [eV] 


tpp [eV] 


NCCO (This work) 
PCCO (This work) 
LSCO (This work) 


1.61 
1.65 
2.76 


8 
8 
8 


1.16 
1.17 
1.41 


0.54 
0.54 
0.66 


NCCO (Ref. [46]) 
LSCO (Ref. |46|) 
LSCO (Ref. [9]) 
LSCO (Ref. [To]) 


1.42 

0.918 

3.5 

3.6 


10 

10 

7.9 

10.5 


1.18 

1.357 

1.5 

1.3 


0.69 

0.841 

0.6 

0.65 



considers the d^z2-r2^ Pz and p-z orbitals, was also con- 
sidered for LSCO. The 6-band Hamiltonian is: 
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(3) 



The hopping parameters tij were obtained by 
downfolding the LDA band structure to six or- 
bitals {dx2-y243z2-r2,Px,y,Pz-z)' The doublc 

counting for the 6-band model is defined as 
Edc = ^d (^d.2-y2 +^ci3.2-r2 - 0-5). Thc samc on- 
site repulsion Ud was considered for the dx2-y2 and 
d3z2-r2 orbitals. 

We solve these models using Dynamical Mean Field 
Theory (DMFT), in which the Green's function is given 
by: [83]. 

Gk(^CJn) = {i^n + /i - Hk - ^{iuJn)y^ , (4) 

where Hk is the Fourier transform of the Ht in Eq. (IT]) 
and Eq. ^. 5] is the self-energy matrix being nonzero 
only in the d orbital. The self energy in Eq. Q is ob- 
tained by solving an Anderson impurity model subject 
to the DMFT self-consistency condition: 



Nk 



kGBZ 



(5) 



where the sum runs over the first Brihouin Zone (BZ), 
and P is projecting the averaged green function onto the 
impurity cluster subspace. 

In this work we use the continuous time quantum 
Monte Carlo impurity solver algorithm 071 |48]. Real 
frequency resolved quantities were obtained by analytic 
continuation of the observables obtained in matsubara 
frequencies. We have crossed checked the analytic con- 
tinuation of the observables obtained on the matsub- 
ara axis using several other numerical solvers: the ex- 
act diagonalization solver (ED) [22 , the density matrix 
renormalization group solver (DMRG) [49 , and the NCA 
solver (Non-crossing diagram approximation). CTQMC 
and ED/DMRG/NCA are complementary tools, working 
respectively on the matsubara and real axis. 

Magnetism was considered within the single site 
DMFT by solving two independent impurity problems, 
while in the case of the 2-site cluster DMFT (c-DMFT), 
the 2-site magnetic unitcell is mapped to a 2-impurity 
cluster. Cluster DMFT improves single site DMFT by 
adding the non-local self energy, not present in the sin- 
gle site DMFT. The 6-band calculation maps the 2- 
correlated orbitals dx2-y2 and d2,z2-r2 to an impurity 
containing two different orbitals. 

Finally, model Hamiltonian ([l]) and ([2| were studied 
previously (for examples see Refs. [27 and [50 ). When 
Ud is large, there is a metal to charge transfer insulator 
transition at integer filling, as a function of the charge 
transfer energy e^ — Cp. However the electron doped 
cuprates fall in a regime where the phases with magnetic 
long range order occupy a large fraction of phase space 
and this regime was not investigated previously. 



PHASE DIAGRAM OF LSCO AND NCCO 

In this section we discuss the phase diagram of NCCO 
and LSCO, and in particular their magnetic properties. 
We compare the various treatments of the short-range 
and long-range correlations. In particular, we consider 
the ordered state within the single site DMFT and within 
the 2-site cluster DMFT. 

In Fig. [l^ we show the magnetic moments of NCCO 
obtained within a three-band description (left panel of 
Fig. [l]). Single site DMFT data of NCCO (blue cir- 
cles, left panel) are in remarkable agreement with exper- 
imental data for ah dopings [51 , though LDA+DMFT 
slightly overestimates the range of stabilization of the 
Neel state: in experiment, the magnetic order of NCCO 
vanishes at doping 5 ~ 0.15, while in our calculations we 
find a slightly larger critical doping of S ^ 0.2. We also 
show data for the cluster cellular DMFT (LDA+cDMFT) 
(grey diamonds). For NCCO, the c-DMFT and the sin- 
gle site DMFT data are almost identical. This is a very 
strong test that the magnetic correlations of NCCO are 
well captured by single site DMFT. In the case of LSCO 



(right panel of Fig. [T]), we found significant difference in 
the region of stability of the magnetic state between the 
single site and cluster DMFT {S < 10%), hence the dy- 
namical short range correlations - absent in single site 
approach - are very important for LSCO in the under- 
doped regime, as reported in [18]. For LSCO, we also 
considered the 6-band theory (the 6-band calculations 
include the description of the p±z and dsz2-r2 orbitals), 
which does not change magnetic moment of the three 
band theory. Experimental values for LSCO [52j are also 
shown and compare qualitatively to our data. 

We also carried out the Hartree Fock calculation, and 
we found that in this static approach the magnetization 
vanishes only at unrealistic large doping S ~ 50% for 
both NCCO and LSCO (dashed line in Fig. [l^), which 
points towards the important role of dynamic correla- 
tions at finite electron and hole doping. Indeed, in the 
Hartree- Fock calculations, the picture is fairly simple: 
The magnetic state is stabilized due to an optimization of 
the Coulomb energy at the expense of the kinetic energy 
(the staggered magnetic order avoids double occupation) 
and in this picture there is no strong difference between 
LSCO and NCCO, as shown in Fig.flk) The same mecha- 
nism is responsible for slight overestimation of the critical 
doping within DMFT. 

When comparing to experiments, it is important to 
keep in mind that two dimensional compounds are not 
able to sustain infinite-range magnetic order at a finite 
temperature (Mermin- Wagner theorem [53]). Therefore, 
the Neel temperature within DMFT should be inter- 
preted as the temperature below which the magnetic cor- 
relations become long but remain finite. This tempera- 
ture can be much higher then the actual Neel tempera- 
ture of the material, which is controlled by the magnetic 
exchange between the two dimensional copper oxide lay- 
ers; and vanishes for a well separated copper oxide planes. 
In Fig. [TJd we show the temperature dependence of the 
magnetic moment for the 3-band calculations of LSCO 
and NCCO. The extracted mean- field Neel temperature 
is about 1500K for NCCO and lOOOK for LSCO, which 
is much higher than the experiment Neel temperature. 
Finally, we emphasize that the magnetic moment of the 
parent compound of LSCO m = 0.42 is larger than the 
one of NCCO m = 0.35, which might suggests that LSCO 
is more correlated than NCCO. The strength of the cor- 
relations can also be quantized by calculating the quasi- 
particle renormalization amplitude Z (see Fig. |2|. First, 
we consider the paramagnetic solutions of both NCCO 
(white circles on the left side) and LSCO (blue circles on 
the right side). We find that the quasi-particle weight is 
going to zero in the parent compound of LSCO: This is a 
signature that LSCO is a charge transfer insulator. For 
NCCO, we find that it is going to a finite value at zero 
doping, which is the signature that NCCO is a paramag- 
netic metal. 

At finite doping, we find that the quasi-particle weight 



of the paramagnetic state is much larger in NCCO than 
in LSCO. This is a signature that NCCO is more metahic 
and less correlated than LSCO. 

We now turn to the calculations for the ordered state 
of both LSCO (red circles on the left) and NCCO (red 
diamonds on the right). Here the average quasi-particle 
weight in the ordered state was obtained by the follow- 



ing formula: Z = ^Pi^^f)/^ 



Pi{^F 



The motivations 



comes from the formula for the specific heat of the mag- 
netic system given by 7^ oc -|^ with i = A^B. 

The quasi-particle weight in the ordered state of 
NCCO is larger than the one for LSCO, which shows 
that the ordered state of NCCO is also less correlated 
than LSCO, and hence the character of the paramagnetic 
state of the parent compounds (paramagnetic metal ver- 
sus paramagnetic insulator) has direct consequences for 
the magnetic solutions. 

We find that the mechanism to open a gap at integer 
filling is totally different for hole and electron doped com- 
pounds, and our results place NCCO and LSCO to dif- 
ferent regions of the ZSA phase diagram. For hole doped 
compound NCCO, the quasi-particles are scattered in- 
creasingly and get a very short lifetime when approaching 
the insulator (charge transfer insulator). In the electron 
doped compound NCCO, the system minimizes its free 
energy by doubling the unit-cell which opens a Slater gap 
(Slater insulator). 

Finally, we note that the presence of magnetism is 
concomitant with an entropy loss , which results in an 
increase of the quasi-particle weight. The direct conse- 
quence is that the specific heat 7 of the ordered state is 
lower than the specific heat of the paramagnet (inset of 
Fig.|2|. 
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FIG. 1: (Colors online) a) We show the theoretical mean- 
value of the staggered magnetization S^ — ^ (n^ ~ "^i) ob- 
tained for LSCO and NCCO by both single site DMFT 
(LDA+DMFT) and cluster cellular DMFT (LDA+cDMFT) 
of the 3-band theory, and DMFT of the 6-band theory. Ex- 
perimental values M{6)/Mo (Mo = M{6 = 0)) for NCCO 
[51] and LSCO [52] are also shown, and for comparison with 
DMFT, we assume Mq = Mdmft{S = 0), where Mdmft 
is the magnetic moment at doping obtained by single site 
DMFT. The general trend compare well to our data, and for 
NCCO we obtain a qualitative agreement for the critical dop- 
ing. For LSCO we find that the magnetization vanish at large 
doping within the single site DMFT approximation, but when 
the possibility of local singlet is allowed (cDMFT) the range 
of stabilization of the long-range order is much reduced. Note 
that the theoretical magnetic moment was obtained at con- 
stant temperature T — 89° K, and experimentally at 2°K for 
LSCO and 8°K for NCCO. b) Magnetization versus temper- 
ature in the parent compounds of NCCO and LSCO obtained 
by single site DMFT for the three band theory. The extracted 
mean-field Neel temperature is about 1500K for NCCO and 
lOOOK for LSCO. 
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FIG. 2: (Colors online) We show the quasi-particle weight 
Z for the 3-band description of LSCO and NCCO. The quasi 
particle weight of the paramagnet (PM) is finite at integer 
filling for NCCO, which is a signature that the paramagnetic 
state of NCCO is a metal. In the ordered state of NCCO 
(AF), the quasi-particle weight is estimated by the specific 
heat of the A and B magnetic sublattices 7^ oc ^^^.^^ and 
i = A, B, and the total quasi-particle weight Z is given by 
Z = Ey^*(eF)/E ^^^- The specific heat 7 of the AF and 

i i 

PM states of NCCO are shown in the inset. 



PHOTOEMISSION AND FERMI SURFACE OF 
NCCO 

The basic quantity describing the electronic structure 
of the material is the electronic spectral function: 



A"(k,cj) 



(ImG(k,a;)"'") 



(6) 



Where a is the orbital index, and k is running through 
the unfolded Brillouin zone. The total spectral weight is 

A(k, cj) = ^ A"(k,cj). Experiments like Angle- Resolved 

a 

Photo-Emission Spectroscopy (ARPES) are able to probe 
the k-dependent spectral functions, and can therefore be 
compared side- by-side with theoretical calculations [84] . 
In this section we investigate the agreement between the 
theoretical spectra and ARPES measurements. 

In Fig. [3^,c,e we show the spectral functions resolved 
in momentum space. In Fig. [3|3,d,f we show the inte- 
grated spectral functions, that show the energy locations 
of the main spectral features. 

In the parent compound (Fig. Isk-b), we find two dis- 
persive peaks separated by the charge transfer gap of 
about 1.2ey, as expected in NCCO [54] . 

The spectral feature below the Fermi level (Fig. Isb) is 
an admixture of oxygens and copper orbitals , commonly 
known as the Zhang- Rice singlet (ZR). It is worth noting 
that our results show that the oxygens orbitals carry no 
magnetic moment. The oxygen sites average the magne- 
tization on its both copper neighbors. 

When the Slater gap opens, there is a spectral weight 
transfer from the upper Hubbard band to the Zhang- Rice 
singlet and the lower Hubbard band (located at — lOeV, 
not shown) , such that minority spectral weight is concen- 
trated in the upper Hubbard band (UHB), and the ma- 
jority spectral weight is mostly present in the lower Hub- 
bard band (LHB) and in the Zhang- Rice singlet (ZR). 

The top of the lower band occurs at (7r/2,7r/2), while 
the bottom of upper band appears at M = (tt, 0), there- 
fore the gap is indirect (see yellow arrow in panel a). 
Those two bands can also be obtained in the simpler 
Hartree Fock approximation, though the size of the gap 
is overestimated in a static mean-field. 

At 10% electron doping (Fig.jsf^-d), NCCO is stih mag- 
netic, and therefore the Zhang- Rice singlet and the upper 
Hubbard band are well separated. Those two features 
are also observed by the simpler Hartree Fock approxi- 
mation. What is clearly not visible in static mean-field, is 
the presence of a very sharp and narrow band slightly be- 
low and above the Fermi level (Fig.lsfc), that corresponds 
to the quasi-particle peak (QP) in the integrated spectra 
(Fig. [SJl). It is worth noting that the optical transitions 
occur at this doping within these narrow bands, from the 
narrow band below the Fermi level to the narrow band 
above the Fermi level, as depicted in Fig. ^ by the ver- 
tical arrow. In the ordered state of NCCO we observe in 
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FIG. 3: (Colors online) Frequency dependent spectral weight 
A(k,w) obtained by LDA+DMFT of NCCO at a) integer- 
filling, c) 10% and e) 20% electron doping for NCCO. A(k,w) 
was obtained along the usual path T — M — K — T in the 
Brillouin zone (see inset of d). The solid lines are the rigid 
LDA bands. The partial density of states of the d and p or- 
bitals are shown in b),d),f). The density of states is showing 
the upper Hubbard band (UHB), the Zhang rice singlet (ZR) 
and the quasi-particle peak (QP) for the doped compounds. 
a) At integer-filling, we find for NCCO that the indirect gap 
is ~ 1.2eV, between (7r,0) and (7r/2,7r/2), as shown by the 
diagonal arrow, and the direct gap is about 1.5eV, as indi- 
cated by the vertical arrow, d) At 10% doping we observe 
a splitting of the quasi-particle peak due to magnetism. The 
splitting of the quasi-particle peak is associated with the mag- 
netic pseudo-gap at (7r/2,7r/2) in panel (c). The horizontal 
arrows pointing from c) to d) are guide to the eyes. The opti- 
cal transitions at 10% doping occur around (tt, 0) within the 
quasi-particle band, as indicated by the vertical arrow, e) 
At 20% doping, magnetism is destroyed and the pseudo-gap 
at (7r/2,7r/2) is closed. The quasi-particle peak in panel (f) 
is clearly related to the spectral weight close to the Fermi 
surface at (tt, 0), as indicated by the horizontal arrow. 



the ordered phase a sphtting of the quasi-particle peak 
(Fig. [sjl) into two structures. The first corresponds to 
the narrow band behow the Fermi level, and its main 
weight is at M = (7r,0), as indicated by the lower hor- 
izontal arrow. The peak slightly above the Fermi level 
is due to the spectral weight at k = (7r/2, 7r/2) (see up- 
per horizontal arrow), and is related to the pseudo-gap 
around k = (7r/2,7r/2). 

Upon larger doping 20% (Fig.|3^-f), magnetism disap- 
pears and the pseudo-gap at k = (7r/2,7r/2) closes. The 
peak slightly above the Fermi level, that was present at 
10% doping, now disappears. 

We now turn to Fig. [4^-d, where we show the spec- 
tral functions resolved in momentum space at fixed en- 
ergies Ef, Ef - O.OSeV, Ef - O.leV and Ef - O.lSeV, 
at finite doping 10%. At this doping NCCO is mag- 
netic. The magnetic Fermi surface in panel (a) has a 
square-like shape structure centered around M = (7r,0). 
At lower energy, in panel d, we observe the presence 
of an arc centered around k = (7r/2,7r/2). This comes 
primarily from the pseudo-gap around momentum point 
k = (7r/2, 7r/2), which is a signature of the magnetic long 
range order (the Fermi surface of the ordered state is 
gapped at /c = (7r/2,7r/2). The Fermi surface of the 
ordered state moves towards the usual Fermi arc shape 
when the system becomes metallic and the pseudogap at 
k = (7r/2,7r/2) is closed in the paramagnet. 

Some aspects of the doped electronic structure can be 
understood in terms of the Hartree Fock rigid band pic- 
ture, for example the holes appear first upon doping at 
the M = (tt, 0) point, but the renormalization of the 
bands, and the multiple peak structure in energy for a 
given momentum point (see Fig. pb) , are not captured in 
static mean-field. 

In Fig. [4^ we show the energy dependence of the spec- 
tral function at a fixed k point (shown in the inset of 
the figure). The peak close to the Fermi energy is con- 
nected to the square-like Fermi surface of panel (a) and 
is hence connected to magnetism. The peak at lower en- 
ergy — 0.2eV, is related to the arc shape of panel d and 
has paramagnetic character. The peak positions are in 
a very good agreement with recent angle resolved photo- 
emission measurements of Ref. [55 also shown in Fig.Efe. 

In Fig.[5^-c we compare side by side experimental data 
(middle panels) [56] and DMFT calculations (upper pan- 
els). The agreement is quantitative showing that our 
approach captures the low energy physics of NCCO. The 
pseudo-gap at (7r/2,7r/2) is also observed in experiments 
[55] (Fig. [sjl, middle panel) in the ordered phase, and 
compares well to our theoretical calculations (Fig. |5]i, 
upper panel). In Fig. [5^ we show the spectral weight 
along the diagonal cut of the Brillouin zone. We observe 
the presence of a sharp kink in the dispersion (waterfall) 
that was also recently reported in experiments [57] . 

In Fig. |6] we show the doping evolution of the Fermi 
surface obtained by theoretical calculations (lower pan- 




FIG. 4: (Colors online) Fermi surface maps obtained by 
DMFT calculations in the ordered state at 10% doping, a) 
at the Fermi energy and b)-d) at lower energies ranging from 
-0.05eV to -0.15eV. The Fermi surface map at the Fermi level 
is shaped by the presence of magnetism, whereas the arc in 
the energy map at energy -0.1 5eV is also present in the para- 
magnetic calculations, e) comparison of A{uj) for a fixed k 
point k = (37r/4, 7r/4) (shown in the inset) obtained theoreti- 
cally (lower curve) and experimentally from Ref. 55 (upper 
curve). The peak at -0.03eV is due to magnetism, and the 
peak at lower energy -0.2eV is associated to the arcs seen in 

d). 



els). We find that upon doping, the Fermi surface moves 
from the square-like magnetic Fermi-surface towards the 
Fermi arcs of the paramagnetic Fermi surface (panels d- 
g). This is explained by the closing of the pseudo-gap at 
k = (7r/2, 7r/2). The agreement with experiments (panels 
a-c) ^58^ is very satisfactory. 
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FIG. 5: (Colors online) Side by side comparison of A{k,uj), 
along the path as depicted in the inset, obtained theoreti- 
cally (upper row) and experimentally (lower row), a)-c) at 
15% electron doping from Ref. ^56], and d) at 13% doping 
from Ref. 55 . e) Comparison between DMFT (left side) and 
experiments of Ref. [57] (right side), along the nodal cut of 
the Brillouin Zone at 17% electron doping. The agreement 
between DMFT and experiments is remarkable. The vertical 
dashed lines are guide to the eyes to illustrate the presence of 
a sharp kink in the dispersion (waterfall) . 



PHOTOEMISSION OF LSCO 

In this section we focus on the spectral functions of 
LSCO. In particular, we focus on the 6-band theory, 
which includes dsz2-r2 and p±z orbitals, which are ex- 
pected to play a role in LSCO, due to the presence of 
apical oxygens. Note that the apical oxygens are absent 
in NCCO. 

In Fig. [7^ we show the momentum resolved spectral 
function of LSCO obtained by DMFT. We observe a di- 
rect gap ~ 1.8eV, which is larger than the gap in NCCO, 
showing that LSCO is more correlated than NCCO. The 
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FIG. 6: (Colors online) Side by side comparison of the Fermi 
surface obtained experimentally (reproduced from Ref. [58] ) 
and obtained by single site DMFT calculations for NCCO in 
the ordered state, a) Experimental results for doping 4%, b) 
for 10% electron doping and c) 15% electron doping. The re- 
sults are compared to DMFT calculations for similar dopings 
d)-g). The Fermi surface at low doping (a) is centered around 
M = (7r,0), and is moving towards the Fermi arc shape (g) 
when magnetism is destroyed. 



partial density of states (Fig. [Tb) shows two dispersive 
features, the upper Hubbard band (UHB) and a the band 
below the Fermi level, the Zhang- Rice singlet (ZR). The 
latter is an admixture of oxygen and copper characters. 
The Zhang-Rice singlet is more incoherent in LSCO than 
in NCCO (see Fig. 1^). 

Fig. [7}3 is a blow up of Fig. [T]^, that displays the in- 
tegrated spectrum on a larger energy scale. The lower 
Hubbard band (LHB) is separated from the upper Hub- 
bard band (UHB) by an energy scale of the order of Ud- 
It is worth noting that the p^^ and dsz2-r2 orbitals have 
a strong weight between — 4eV and — leV, and the in- 
plane oxygens are located at — 5eV. Hence, the additional 
orbitals p±z and dsz2-r2 hybridize with the Zhang-Rice 
singlet, and change the theoretical description of LSCO 
for energies larger than leV. 

In Fig. [8^ and [8]:: we focus on the momentum resolved 
spectral functions of doped LSCO. Figures ^jp and [8]d 
are the corresponding integrated spectral functions. At 
10% doping (Fig. Isb), the Zhang- Rice singlet has an in- 
coherent contribution (ZR) and a coherent part - the 
quasi-particle peak (QP). The coherent part (QP) is the 
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FIG. 7: (Colors online) a) Frequency dependent spectral 
weight A(k,w) obtained by LDA+DMFT of a 6-band model 
description of the parent compound of LSCO. b) Partial den- 
sity of states of the dx2-y2 , d3z2-r2 , P(x,2/) and p±z orbitals. 
We observe a direct gap of 1.8eV in LSCO. Notice that the 
spectral weight is very incoherent close to the Fermi energy 
in the lower band, c) Partial density of states on a larger en- 
ergy scales. The lower Hubbard band (LHB) is located at a 
very low energy — lOeV, and the upper Hubbard band (UHB) 
is also shown. The d3z2-r2 and p±z orbitals have a strong 
weight between — 4eV and — leV. 



around a nodal cut of the Brillouin zone. We observe the 
presence of a very sharp kink in the dispersion (water- 
fall)^ in agreement with experimental data of Ref. [59] 
reproduced in Fig. |8ji. 

Fig. [9^-d displays the orbital resolved spectral func- 
tions, obtained within a 6-band theory, for various dop- 
ings on a wide energy scale. Upon doping (Fig. [9|3-d) we 
observe the presence of the quasi-particle peak, and the 
upper Hubbard band (UHB) smears out. Additionally, in 
the doped compound the dsz2-r2 and p±z orbitals have 
a larger weight close to the Fermi level. 

Fig. [9^ show the relative occupation (in hole notation) 
of the dsz2-r2 and p±z orbitals upon doping. Our results 
are in agreement with Ref. [60 and we find that even 
after inclusion of the apical oxygens, there is no satu- 
ration of the occupancy observed around doping 0.2 in 
the x-ray absorption spectroscopy (XAS) experiments of 
Ref. [61 . We notice however that LDA+DMFT does 
capture the evolution of the ratio of the occupancies of 
apical and planar oxygens: For doping S < 20%, the 
holes go mainly to the dx2-y2 and P{x,y) orbitals, and for 
larger doping S > 20%, the holes start to fill the dsz2-r2 
and p±z . In Fig. [9]F we report the experimental data of 
Ref. [62] for side-by-side comparison with DMFT calcu- 
lations. Our calculations reproduce the rapid increase of 
the occupancy of Pz around doping 0.2. 

We note that modeling XAS more accurately may re- 
quire downfolding the LDA bands over a much larger 
energy range to include more orbitals, or the inclusion of 
the doping dependence modifications of the apical oxy- 
gen distance to the copper oxide layer pointed out in Ref 
[63 . A more accurate modeling of the XAS, including 
the core hole potential, as it was done for the core level 
photoemission in Ref [64], might also be necessary. 



narrow band below the Fermi level in Fig. [8^. The ver- 
tical arrows in Fig. [8^ and [Sj^ highlight the location of 
the direct transitions from occupied states to unoccupied 
states. These transitions are important for the optical 
conductivity (which we discuss in the next section). 

For comparison, we also show the theoretical descrip- 
tion of LSCO without the apical oxygens in Fig. [8^-f. 
The main difference between the 3-band and the 6-band 
descriptions, is that the incoherent part of ZR is narrower 
in the 3-band theory. The vertical transitions, marked 
with vertical arrows in Fig. |8]3,e, highlight large contri- 
butions to the optical conductivity, and one can see that 
the optical transitions occur at different energies in the 
two models. 

Fig. [Sb is the momentum resolved spectral function 
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FIG. 9: (Colors online) Spectral functions for the 6-band the- 
ory of LSCO in a) the parent compound and at b) 10% and c) 
20% doping and d) 30% doping. The dx2-y2 , d3z2-r2 , P(x,y) 
and p±z components are shown. Note that the lower Hubbard 
band (LHB) at — 9eV is both of dx2-y2 and d3z2-r2 charac- 
ter. At finite doping (b) and (c), there is a quasi-particle peak 
(QP) close to the Fermi energy. Note that the d3z2-r2 and 
p^2 weight gets closer to the Fermi level for higher doping (c). 
These latter orbitals will therefore be important to describe 
high doping and/or high energy spectroscopy, e) Ratio of the 
hole occupancy of the d3z2-r2 and dx2-y2 orbitals (squares) 
and of the p±2 and V{x,y) orbitals (circles) obtained by DMFT, 
and f) experimental data are also shown for comparison _62] . 
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FIG. 8: (Colors online) a),c) Frequency dependent spectral 
weight A(k,w) obtained by LDA+DMFT of a 6-band model 
description of doped LSCO. b) and d) Partial density of states 
of the dx2-y2 , d3z2-r2 , V{x,y) and p±;z orbitals. a) and b) are 
obtained in the ordered state, e) and f) are obtained within 
a 3-band description of doped LSCO and are shown for com- 
parison. The optical transitions occur at different energies for 
the 3-band and 6-band theories (see vertical arrows in c and 
e). Also note that the weight of the d3z2-r2 and y>±z orbitals 
in b) and d) close to the Fermi level is consequent, which jus- 
tifies that a 6-band description is necessary. The solid lines 
in a),c) and e) are the rigid LDA bands, g) A(k,w) along the 
nodal cut of the Brillouin Zone for the paramagnetic 6-band 
theory at 10% doping, compared side by side with h) exper- 
imental data of Ref. 59 . The dashed line are guide to the 
eyes to illustrate the sharp kink [waterfall) in the dispersion. 



13 



OPTICAL PROPERTIES OF LSCO AND NCCO 

We now turn to the optical conductivity. It was previ- 
ously computed for the 3-band model of LSCO [31] and 
we now generalize the results for the 6-band description 



of LSCO and we also compute it for NCCO (see Fig. 10). 
The optical conductivity in LDA+DMFT is given by: 



^ ^ Nk^ he J 



f{x-uj) - f{x) 



X Tr (pka {X - CJ) Vk/3ka (^) Vk) (7) 

Where c is the interlayer distance, and the density matrix 
p is defined by 



/5k.(x) = ^(GL(a;)-Gk.(a;)) 



(8) 



The bare vertex for a multiple orbital problem v? = 



dm 



dkx 



^ {Qx ~ Qx) H^ is obtained following the steps 
of Ref [65]. The Peierl phase i {q^ — Qx) H^ plays an 
important role in particular for the ordered state (as dis- 
cussed in Ref ^65j, if this phase is not considered the 
optical conductivity depends on any artificial folding of 
the Brillouin zone). 



In Fig. 10 i we show the theoretical optical conduc- 
tivity of NCCO at integer filling (red curve) and at 10 
percent doping (blue curve). The undoped compound 
has a sharp onset at an energy of the order of 1.5 eV 
which we interpret as the direct gap (slightly larger than 
the charge transfer gap in Fig. ^jp. 

Doping introduces several new features (blue line in 
Fig. flOk). The l.beV optical peak disappears and the 
weight is transferred to lower energy in the form of a 
Drude peak and a mid infrared peak at cj ^ 0.2eF. The 
optical conductivity also displays a peak in the magnetic 
solution at a much lower frequency u ^ 0.035eV (see left 
inset of Fig. 10 i). 

Below 0.5 eV, vertical transitions are among the quasi- 
particle bands of the magnetic DMFT band structure. 
This involves a continuum of k points, but it is likely to 
be controlled by saddle points in the reciprocal space. 
One saddle is at M = (tt, 0) and transitions close to that 
point, indicated by a vertical arrow in Fig. ^, give rise 
to the peak in the optical conductivity at 0.2 eV. Tran- 
sitions close to the midpoint between K = (tt, tt) and 
r = (0, 0), give rise to the small peak at 0.035 eV in the 
optical conductivity. Both peaks are characteristics of 
the quasi-particle band structure in the magnetic state of 
NCCO, and these features are not present in the absence 
of magnetic order, as shown in Fig. [TT^ . The absence of 
the peak at 0.035eV and 0.2eV in the optical conductiv- 
ity is explained by the absence of vertical transition in 
the k dependent spectral weight A{k^uj) (see Fig. [lib), 
which are present in the ordered state and highlighted by 
the vertical arrow of Fig. ^jp. 



The quasi-particle peak in the ordered state is split 
into two narrow bands, one above and one below the 
Fermi level (see Fig. [3J1), while the splitting is absent 
in the paramagnetic state (see Fig. fllfc) . This is due to 
the absence of the pseudo-gap around (7r/2,7r/2) in the 
paramagnetic state (see Fig. [lip). As a consequence, 
the paramagnetic Drude peak is featureless, while in the 
ordered state we observe several low energy peaks in the 
optical conductivity (at 0.035 eV and 0.2 eV). 

The agreement between DMFT and experimental data 



J66[ [67] (dashed line of Fig. 10 i) is qualitative and our 
theory connects the peak in the experimental optical con- 
ductivity at 0.4 eV with magnetism, in agreement with 
Ref. [68]. Additionally our study allow us to connect 
this peak with the spectral weight below the Fermi en- 
ergy at the M = (tt, 0) point, which is present in both 
the paramagnetic and the ordered state. 

We note finally that the peak in the optical conductiv- 
ity at smaller energy 0.035eV is observed at this energy 
in experiments (see Fig. 2 of Ref. [69 ), which is only 
present within the ordered state and disappears at higher 
temperature in the paramagnet, which suggest that this 
feature is not connected to phonons. Our calculation con- 
nects this peak to the presence of magnetism and to the 
pseudo-gap around (7r/2,7r/2) (see Fig. [3]::). 

We now turn to the description of the optical conduc- 
tivity of LSCO. In Fig. [TOJd we show the optical con- 
ductivity of the 6-band description of LSCO in the par- 
ent compound (red curve), at 4% (black curve) and 16% 
doping (blue curve). For comparison we also show ex- 
perimental data of Ref. [70 at 0% doping (dotted line), 
6% (short dashed) and 20% (long dashed). Note that 
there is a quantitative agreement between the theory and 
experiments for frequencies smaller uj < 2eV. At larger 
frequencies u > 2eV, the optical conductivity in experi- 
ments is larger than the theoretical one, which is due to 
optical transition from additional orbitals that are not 
present within our calculations, and which contributes 
at high energies. 

For comparison, we now also consider the three band 



theory of LSCO. In Fig.lTO):; we show the optical conduc- 
tivity of the magnetic state of the parent compound (red 
curve) and of the doped LSCO at 20% doping (yellow 
curve). We also show the experimental data of Ref. [70] 
at 0% doping (dotted line) and 20% doping (dashed line). 
The vertical arrow emphasizes the disagreement between 
theory and experiments. The strong differences between 
the 3-band and 6-band calculations for the optical con- 
ductivity of the parent compound is related to the strong 
weight of the p±z and dsz2-r2 orbitals below the Fermi 
level, which are mainly occupied and contribute signifi- 
cantly to the optical conductivity. This weight is obvi- 
ously absent in the three band theory. Note also that 
at finite doping there is a peak at 0.8eV in the opti- 
cal conductivity of the 3-band theory, which is absent in 
experimental data. This peak is related to the optical 
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transitions shown in Fig. |8^ by the vertical arrows. This 
artifact of the 3-band calculations is cured by introducing 
the dsz2-r2 and p±z orbitals. 

To quantify the rate of the redistribution of optical 
spectral weight, we computed the effective electron num- 
ber per Cu atom defined by 



^i// 



2meV 



a{(jo)d(jO^ 



(9) 



where rrie is the free electron mass, and V is the cell vol- 
ume containing one formula unit. N^ff is proportional 
to the number of electrons involved in the optical excita- 
tions up to the cutoff A. 

Our results for N^ff are displayed in Fig. p!2]and com- 
pared to experimental data taken from Ref. [71], [72] and 
[73] . Notice a favorable agreement between the theory 
and experiment, for which the use of the realistic elec- 
tronic structure is essential. Moreover, we also show in 



Fig. 12 the static mean-field results. We emphasize that 
the static mean field is not able to reproduce qualitatively 
the experimental trend. This highlights the importance 
of the dynamical fiuctuations not taken into account in 
Hartree Fock. The trend of N^ff is qualitatively simi- 
lar for both LSCO and NCCO, and therefore does not 
show that LSCO and NCCO are qualitatively different: 
the former being a Mott and the latter being a Slater 
insulator. 

We also compute the temperature dependence of A^e// 
for NCCO and LSCO at a fixed density. In Fig. [T3| l we 
show the theoretical results for NCCO. We plot in the 
same figure the temperature dependence of A^e// (left 
scale) and the temperature dependence of the magnetic 
moment (right scale). We find that N^ff is reaching a 
maximum value when the magnetization is suppressed 
by thermal fiuctuations at ~ 400K. There is actually a 
variation of Neff inside the ordered phase, that can be 
explained by the destruction of the magnetic pseudo-gap. 
Once the magnetization is quenched by the temperature, 
heating the system further reduces Neff. The qualita- 
tive agreement with experimental data extracted from 
Ref. [72 is very encouraging. 

The temperature dependence of Neff for LSCO is 
shown in Fig. [13)3. We emphasize that the temperature 
dependence of Neff shows an opposite trend for LSCO. 
When doping the parent compound, there is a decrease 
of Neff in LSCO, whereas there is an increase of Neff in 
NCCO. Therefore the temperature dependence of Neff 
clearly shows a distinct behavior for a Slater and a Mott 
insulator. Note that the same general trend of the tem- 
perature dependence of Neff is observed in experiments 
for LSCO, but for a larger doping than the one consid- 
ered here (optical data for 13% doping in Fig. 6 of Ref. 

124,)- 

The trend of Neff can be understood in a simple pic- 
ture. In a Slater picture, the onset of antiferromagnetism 
reduces the Coulomb correlations (double occupancy) at 
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FIG. 10: (Colors online) a) Theoretical optical conductiv- 
ity of NCCO in Q~^cm~^ at integer- filling (red line). We 
find that NCCO has an optical gap of about 1.5eV, which is 
larger than the direct gap ^ 1.2eV. We also show data for 
doped NCCO (blue line). Note that for doping smaller than 
20%, we observe a small peak at a small energy ^ O.OSSeV 
(see the inset) . We also observe a second peak at larger energy 
scale ~ 0.2ey , which corresponds to optical transitions within 
the Zhang-Rice singlet (see Fig. [g]^). For comparison we also 
show the infrared optics of Refs [66l [67] (dashed black line) . 
b) Optical conductivity of the 6-band theory of LSCO for the 
parent compound (red line) and doped LSCO. For compar- 
ison we also show experimental data of Ref. [70 (short and 
long dashed lines) . c) Optical conductivity of the 3-band the- 
ory of LSCO in the parent compound (red line) and doped 
LSCO (yellow line), and experimental data of Ref. [70 (short 
and long dashed lines). The vertical arrow in panel (c) em- 
phasizes the disagreement between theory and experiments. 
Note also that at finite doping there is a peak at 0.8eV in 
the optical conductivity of the 3-band theory, which is absent 
from experimental data. This peak is related to the optical 
transitions shown in Fig. ^b (see the vertical arrows). All 
calculations were done in the ordered state. 
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FIG. 11: (Colors online) a) Theoretical optical conductiv- 
ity of paramagnetic NCCO at 10% doping is shown. In the 
paramagnetic state of NCCO, there are no peaks at 0.035eV 
and 0.2eV which are obs erve d in the magnetic state of NCCO 
at same doping (see Fig 10 a), b) We show the k dependent 
spectral weight A{k, cj) for NCCO at 10% doping, c) We show 
the corresponding total density of states. In the paramagnetic 
state of NCCO there is no splitting of the quasi-particle peak 
QP (seeFig.lsli). 



the expense of the kinetic energy. The opposite is true 
in a charge transfer insulator. Consequently, in a Slater 
insulator the kinetic energy becomes less negative as the 
temperature decreases while the opposite happens in a 
charge transfer insulator. The kinetic energy as a func- 
tion of temperature is readily available in the theory and 
is displayed in Fig. [14^ and Fig. \i^ : and the Coulomb 
energy is displayed in Fig. [l4}3 and Fig. [l4]d. 

In Fig. [T4]F we display the imaginary part of the self 
energy at zero frequency, ImT^iuj = 0), as a function of 
the temperature, which is proportional to the scattering 
rate A = —ImT^{u = 0). We find that the scattering 
rate of LSCO is strongly reduced for T < T^eei {TNeei is 
highlighted by the red vertical arrow), which is consistent 
with the kinetic energy reduction observed for LSCO at 

The scattering rate of 
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FIG. 12: (Colors online) We show the dimensionless in- 
tegrated optical conductivity Neff for the 3-band study of 
LDA+DMFT done on LSCO 31 and NCCO, obtained with 
a cutoff A = 1.2 (A = 1.5) for NCCO (LSCO). Experimen- 
tal data for LSCO (red circles 7X and yellow squares 73 ) . 
and NCCO [72] (open diamonds) are shown. The increase of 
Neff is similar for both compounds. The dashed line indicates 
the results obtained by the Hartree Fock approximation. We 
also show for comparison results obtained within the 6-band 
theory (green squares). The agreement between theory and 
experiments is quantitative. 



NCCO is much smaller and is only weakly temperature 
dependent, which is a consequence of NCCO being less 
correlated than LSCO. 

Hence the location of NCCO and LSCO, relative to 
the metal-to-charge-transfer-insulator boundary |TT], ac- 
counts for the observed trends in the temperature depen- 
dence of the optical conductivity. Similar trends of tem- 
perature dependence of the kinetic energy for both elec- 
tron and hole doped cuprates were reported in Ref. [75]. 



Fig. 15 i displays the Kinetic energy difference between 
the paramagnet and the ordered state, as a function of 
doping at a fixed temperature T=89K, for NCCO and 
LSCO. The difference of Coulomb energy is shown in 
Fig. [15)3. We find that in NCCO, the ordered phase is 
stabilized by optimizing the on-site repulsion (Coulomb 
energy), at a cost in the kinetic energy. This is typical 
for a Slater insulator, and is also captured by a simple 
Hartree Fock static mean-field theory. In LSCO the trend 
is opposite, the ordered phase is stabilized by having a 
lower kinetic energy, at a cost in the Coulomb energy. 

Very interestingly, in a one band theory an increase of 
the order parameter (S^) leads necessary to a decrease 
of the Coulomb energy [76 : 



{{Sn') = \{{n,-n,f) = \{n-2{n^n^)) 



And therefore 



low temperature (see Fig. 14 



(ntn;) = |-2((5^f) 



(10) 
(11) 
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FIG. 13: (Colors online) a) 3-band theoretical normalized 
variation of Neff, ANeff = A^e//(T) - Neff(T = 89° K), 
at 10% electron doping for NCCO (red circles, left scale). 
Neff is reaching a maximum value when the magnetization 
(orange squares, right scale) is destroyed by the thermal fluc- 
tuations. The decrease of Neff at low temperature can be 
explained by the opening of a pseudo-gap in the ordered 
phase. The data were obtained by single site DMFT cal- 
culations in the ordered phase. The dashed line corresponds 
to experiments (see Fig. 7c of Ref [72]), where they measured 
Neff (A = O.OSey) (contribution due to the Drude peak) and 
Neff {A = O.SeV)- Neff {A = 0.2ey) (contribution due to the 
pseudogap). The dashed line corresponds to the sum of these 
two contributions, b) 3-band single site DMFT of the ordered 
phase of LSCO is shown for 5% doping. Note that the trend 
of Neff is opposite between LSCO AND NCCO, which is a 
signature that NCCO is a Slater insulator and LSCO a Mott 
insulator. 



In the ordered phase there is hence an increase of {S^) 
and a decrease of {n^n^) for a fixed density n or doping 
5. 

In a 3-band theory this is not the case, since the den- 



sity n in formula (11) is not the total density but the 
density of the d orbital Ud- The latter quantity is not 
fixed at a given doping, and is hence increased when the 
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FIG. 14: (Colors online) 3-band theoretical temperature de- 
pendence of the a) Kinetic energy 1-Lt (equation [l]) and b) 
Coulomb enery 1-Lu (equation [2]) of NCCO at 10% doping, 
and c) Kinetic and d) Coulomb energy of LSCO at 5% dop- 
ing. The red area highlights the temperature region where 
the solution is magnetic (AF), and the solution is paramag- 
netic (PM) in the blue area. a)-b) are showing that there is 
a kinetic energy optimization when LSCO becomes an anti- 
ferromagnet, which is proper to the Mott insulator, and c)-d) 
show that NCCO is a typical Slater insulator, which opti- 
mizes the Coulomb (local onsite repulsion) energy when it 
becomes an antiferromagnet, at the expense of a worse Ki- 
netic energy. This is consistent with the theoretical Neff 
shown in Fig. 13 f) We show the imaginary part of the self 
energy at zero frequency ImT^{uj — 0) in function of the tem- 
perature. The scattering rate of LSCO A = —ImT^{uj = 0) 
(red lines) is strongly reduced for T < TncbI (TNeei is high- 
lighted by the red vertical arrow), which is consistent with 
the kinetic energy reduction observed for LSCO at low tem- 
NCCO (blue lines) is showing a small 



perature (Fig. 14 



scattering rate weakly dependent on the temperature, which 
is consistent with NCCO being more metallic than LSCO. All 
calculations were obtained by CTQMC in the ordered state. 
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FIG. 15: (Colors online) 3-band theoretical energy differ- 
ences between the antiferromagnetic and the paramagnetic 
phases. We show the doping dependence of the a) Kinetic en- 
ergy Ht (equation [l]) and b) Coulomb energy Hu (equation [2]) 
of NCCO and LSCO at fixed temperature T=89K. There is a 
kinetic energy optimization when LSCO becomes an antifer- 
romagnet, which is proper to the Mott insulator, and NCCO 
is a typical Slater insulator, which optimizes the Coulomb 
(local onsite repulsion) energy when it becomes an antifer- 
romagnet, at the expense of a worse Kinetic energy. This is 
consistent with the temperature dependence of the theoretical 
Neff shown in Fig. |13| and with the temperature dependence 
of the Kinetic and Coulomb energies |14| All calculations were 
obtained by CTQMC. 



magnetization is increased. For the parent compound of 
LSCO, we found that magnetic correlations lead to an 
increase of Ud of 1.4%, and to an increase of the double 
occupancy Ud^ridi of 6%. Hence our results highlight a 
significant difference between the single band and 3-band 
theoretical description of LSCO and NCCO. 



STRENGTH OF CORRELATIONS IN LSCO AND 
NCCO 

We finally extended the 3-band calculations to other 
values of e^^ — e^, in order to study the dependence of our 
results on the charge transfer energy Cd — Cp. We em- 
phasize that the charge transfer energy plays the role of 
an effective repulsion U in the one band model language, 
and hence controls the strength of the correlations in a 
3-band theory. For instance, it was shown in the seminal 
ZSA paper [11] that if the Coulomb repulsion of the d 
orbital is larger than the charger transfer energy e^^ — e^, 
the size of the gap in the paramagnet is independent of 
the Coulomb repulsion, and the strength of correlations 
is set by the charge transfer energy. 

In order to study the strength of the correlations, we 
computed the jump in the chemical potential (5/i for in- 
finitesimal doping of the parent compound. This quan- 
tity gives an estimation for the gap in the spectral func- 
tions of the parent compound, which is around 1.2 eV 
and 1.8 eV in NCCO and LSCO respectively. We present 
the data for both the ordered state and the paramagnetic 
state of LSCO and NCCO in Figs. [16^) and [iT^), respec- 
tively. The jump of chemical potential in the paramagnet 
gives an estimation of the critical charge transfer energy 
Ac2, which is the maximum charge transfer energy that 
allows the metallic solution. However at Ac2 the ordered 
solution has a substantial gap, which is closely related to 
6 /J. in the ordered state. We find that S/j. in the ordered 
state of NCCO is around 1.2 eV and 1.8 eV in LSCO. 



Hence, we conclude that NCCO (16 a) is slightly below 
the charge-transfer-insulator-to- metal transition Ac2- In- 
deed, below Ac2, the magnetic long-range correlations are 
necessary to open a gap, and slightly above Ac2 the para- 
magnetic gap is much smaller than the slater gap induced 



by the nesting, as shown in Fig. [I6p. 

The two solutions of the DMFT equations for NCCO 
are shown in panels [T6)3,c. The first one is non-magnetic 
and metallic, and describes a material in the absence 
of long range order. The second is insulating and an- 
tiferromagnetically ordered, with a charge transfer gap 
of 1.2eV. Since the non-magnetic solution is metallic, the 
magnetic long range order is responsible for the insulat- 
ing nature of NCCO (Slater insulator). 

We find that LSCO (Fig.[l7^) is above Ac2, as reported 
in a recent study [31]. Indeed, the parent compound of 
LSCO is only weakly affected by the presence of mag- 
netic order, the size of the gap is only slightly increased 
when magnetic order is present. The DMFT equations 
for LSCO have two solutions, shown in panels [17)3, c. The 
first one is paramagnetic and the second is antiferromag- 
netically ordered, with a charge transfer gap of 1.8eV. 
Since the paramagnetic solution is insulating, the mag- 
netic long range order is not responsible for the insulating 
nature of LSCO (Mott insulator). 



In Fig. 18 i-c, we show the density of states of the 3- 
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FIG. 16: (Colors online) a) We show the jump in the chem- 
ical potential b\i in the ordered state of NCCO (blue circles) 
for other values of the charge transfer energy e^ — e^. The 
results are obtained for the 3-band single site DMFT. For 
comparison, we also show the jump in the chemical potential 
of the paramagnet (red squares). There is a quantum criti- 
cal point Ac2 for the paramagnetic state, with respect to the 
charge transfer energy, corresponding to the metal to charge 
transfer insulator transition. The physical value obtained by 
LDA+DMFT for td - Cp places NCCO below the boundary, 
b) Density of states of the paramagnetic state of NCCO. The 
quasi-particle peak close to the Fermi energy (QP) is a signa- 
ture that NCCO is paramagnetic metal, c) Density of states 
of the ordered state of NCCO. 




FIG. 17: (Colors online) a) We show the jump in the chemi- 
cal potential S/n in the ordered state of LSCO (blue circles) for 
other values of charge transfer energy Cd — Cp. The results are 
obtained for the 3-band single site DMFT. For comparison, 
we also show the jump in the chemical potential of the para- 
magnet (red squares). There is a quantum critical point Ac2 
for the paramagnetic state, with respect to the charge transfer 
energy, corresponding to the metal to charge transfer insula- 
tor transition. The physical value obtained by LDA+DMFT 
for Cd — Cp place LSCO above the boundary, b) Density of 
states of the paramagnetic state of LSCO, which shows the 
presence of a gap. This is a signature that LSCO is a Mott 
insulator c) Density of states of the ordered state of LSCO. 
The gap in the density of states is of similar size for the para- 
magnetic insulator and the ordered state as shown in details 
in Fig. fTsl 
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FIG. 18: (Colors online) Density of states of the 3-band de- 
scription of the parent compound of LSCO in the paramag- 
netic (dashed line) and ordered state (red area) obtained with 
different solvers : a) exact diagonalization (ED), b) CTQMC 
and c) density matrix renormalization group (DMRG). All 
data show that LSCO is a paramagnetic insulator, and that 
the size of the gap obtained by the density of states is sim- 
ilar for both the paramagnet and the ordered state (within 
~ 10%). d) We show the variation of the doping 6 with re- 
spect to the chemical potential /i for both the ordered state 
(SDW) and the paramagnetic state (PM) obtained by exact 
diagonalization (ED) and CTQMC. There is a jump in the 
chemical potential 6/j. (plateau at ^ = 0) of similar size for 
all data. This shows that the magnetic correlations do not 
strongly affect the insulating properties of LSCO. 



band description of the parent compound of LSCO using 
various numerical tools. In panel a we show the density 
of states obtained by ED for the paramagnet and the or- 
dered states. In panel b we show the density of states 
obtained by CTQMC and in panel c the density of states 
obtained by a recent DMRG solver [49]. We conclude 
that the gap in the ordered state of LSCO is of similar 
size than the gap obtained in the paramagnet, indepen- 
dently of the numerical solver used to solve the DMFT 
equations. In Fig. [iSJl we show the doping as a function 
of the chemical potential, and it displays a plateau re- 
lated to the charge gap at integer filling. The agreement 
between the different solvers gives us confidence in these 
results. 

The asymmetry between both NCCO and LSCO, being 
below and above the charge-transfer-insulator-to-metal 
transition Ac2, is a simple explanation for the strik- 
ing asymmetry in their phase diagram. For LSCO, the 
magnetic correlations are destroyed rapidly upon doping, 
while in the NCCO they survive up to high doping, as 
shown in Fig. [l^. 

In this section, we emphasized that the physical origin 
of the asymmetry between LSCO and NCCO lies not 
only in the different values of the oxygen-oxygen overlap, 
which controls the curvature of the Fermi surface, an 
effect that is captured in model Hamiltonian studies, but 
also in the different values of the charge transfer gap in 
the two structures. The latter has an electrostatic origin, 
the electron doped material lacks the negatively charged 
apical oxygen, which increase the electrostatic potential 
at the copper site. 

For comparison, we now also discuss results of the one 
band Hubbard model. In Fig. [19] we show the gap in the 
density of states for the paramagnetic insulator and for 
the antiferromagnetic insulator. In the one band model 
exists a critical point Uc2 which separates the paramag- 
netic metal at small U/t from the paramagnetic insulator 
at large U/t (squares). The magnetic solution is always 
insulating (circles). Using typical values for the hopping 
parameter t [77](t = 0.42eV for NCCO and t = 0.43eV 
for LSCO), and typical values for the gap A in the or- 
dered state for LSCO and NCCO (A = 1.2eV for NCCO 
and A = 1.8eV for LSCO) we can locate the compounds 
in the one-band model phase diagram. Both LSCO and 
NCCO are below Uc2 in this picture, in agreement with 
Ref [30 . There is therefore a strong difference in the 
physical conclusions obtained by the one-band calcula- 
tions and the LDA+DMFT. 



CONCLUSION 

In conclusion, we carried out a comparative study of 
NCCO and LSCO using a realistic LDA+DMFT ap- 
proach. The LDA+DMFT studies achieved overall good 
agreement with optical conductivity and ARPES studies 
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FIG. 19: (Colors online) We show calculations done on the 
one band Hubbard model for various Coulomb repulsion U 
and transfer integral t ratios. In the one band theory there is 
also a quantum critical point Uc2 that is the minimal repulsion 
driving an paramagnet to an insulator. Locating the com- 
pounds LSCO and NCCO by fitting the gap A (A = 1.2eV 
and t = 0.42eV [77 for NCCO and A = 1.8eV and t = 0.43eV 
for LSCO), we place both LSCO and NCCO below Uc2 in the 
one band picture. It is worth noting that for large enough 
U/t the size of the paramagnetic gap is close to the gap in 
the ordered state. Results were obtained by using single site 
DMFT with exact diagonalization. For comparison we show 
theoretical calculations of Ref. [37 , done for the paramagnet 
(long dashed line) and the ordered state (short dotted line), 
and the theoretical calculations of Ref. 36 in the param- 
agnetic state (white circle) and in the ordered state (black 
circle) . 



in a broad range of dopings and a wide range of energy 
scales up to energies of the order of the charge transfer 
energy. 

In particular, for NCCO we found that static mean- 
field theory is not sufficient to describe many qualitative 
features due to the presence of multiple peaks in the elec- 
tronic spectra. The description of these features require 
more sophisticated methods and studies along those lines 
were carried out for NCCO in Refs [^[[78]. We demon- 
strated that LDA+DMFT successfully describes these ef- 



fects in NCCO. 

While single site DMFT is already a good methodology 
to describe the phase diagram of NCCO, cluster correc- 
tions are important in the underdoped region of LSCO. 
This indicates the importance of singlet formation or the 
possible importance of other ordered states in this region. 
However, our LDA+DMFT gave a remarkable agreement 
with experiments for the optical conductivity, for doping 
smaller than 20% and energy scale uj < 2eV, when the 
<^3z2-r2 and p±z are considered. For larger doping and 
higher frequencies, additional bands must should consid- 
ered for a proper description of the photoemission spec- 
tra. 

Moreover, in agreement with Ref ^60^ we found that 
even after inclusion of the apical oxygens, single site 
DMFT does not capture the saturation of the occupancy 
observed around doping 0.2 in the x-ray absorption spec- 
troscopy (XAS) experiments of Ref. [61 . We notice 
however that LDA+DMFT does capture the evolution 
of the ratio of the occupancies of apical and planar oxy- 
gens, including a rapid increase of the occupancy of Pz 
around doping 0.2. We note that modeling XAS more ac- 
curately may require: i) Taking into account additional 
LDA bands and more orbit als, ii) To include the doping 
dependence of the apical oxygen position pointed out in 
Ref [63|, iii) Non-local correlations beyond DMFT might 
be relevant as suggested in Refs [60l ITOHST] , iv) To in- 
clude the core hole potential, as it was done for the core 
level photoemission in Ref ^64^, might also be necessary. 

We achieved a successful description of many nor- 
mal state physical properties (optical conductivity, 
ARPES, stability of magnetism) of these two archety- 
pal cuprates. The overall quantitative agreement be- 
tween LDA+DMFT and experiments, for both LSCO 
and NCCO gives us increased confidence in LDA+DMFT 
as an approach to strongly correlated materials. 

We emphasize that it did not require the introduction 
of ad hoc parameters, such as a doping dependent in- 
teraction strength. The sensitivity of the results on the 
choice of the charge transfer energy, which essentially de- 
termines the strength of the correlations, was explored. 
We emphasize that the charge transfer energy for NCCO 
and LSCO is obtained from ab initio calculations, and 
for the values obtained by first principle calculations we 
obtained good agreement with experiments. 

It is remarkable to obtain good quantitative agree- 
ment between theory and experiment, on energy scales 
of the order the charge transfer energy. For example, 
our methodology gives a good agreement with experi- 
ments for the integrated spectral weight of the optical 
conductivity with a cutoff above the charge transfer gap 
in LSCO. For this quantitative agreement it is essential 
to use a multiband model including the apical oxygen or- 
bitals. The fact that these are absent in NCCO accounts 
for the good agreement with the magnitude of the opti- 
cal conductivity obtained for this material in our earlier 
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publication [32] , 

Moreover, our methodology applied to two different 
materials captures not only their similarities, as for ex- 
ample the doping dependence of the integrated spectral 
weight, but also their differences, as the doping depen- 
dence of the magnetic moment. These differences stems 
from their different electronic structure and is well cap- 
tured by LDA+DMFT. Differences between the electron 
and hole doped cuprates have been noticed by many au- 
thors. They were interpreted as arising from both the 
bare hopping integrals in the one band model Hamil- 
tonian and the Coulomb interaction U in the Hubbard 
model. 

However, LDA+DMFT and multi-band theories are 
also able to capture a more fundamental difference re- 
sulting from the different strength of the correlations in 
these two materials, which is driven by the charge trans- 
fer energy. Using our methodology, we determined the 
relative strength of the correlations of LSCO and NCCO. 
LSCO was found to be on the insulating side of the ZSA 
phase boundary confirming the preliminary conclusions 
of reference [31] , but in disagreement with Ref [30^^ which 
classifies all the parent compound of the copper oxides as 
Slater insulators. The electron doped compound NCCO 
was found to be on the metallic side of the ZSA phase 
boundary ^ in agreement with ref [30^ , but in disagree- 
ment with Ref [46]. Indeed, the lack of apical oxygens 
in NCCO reduces the charge transfer energy relative to 
LSCO, placing these two materials on two different sides 
of the ZSA phase boundary. 

This work is complementary to our earlier work [32] . 
placing the copper oxides in a region of intermediate cor- 
relation strength. For materials in this region of parame- 
ters, we found that the location of the material relative to 
the ZSA boundary has important physical consequences. 
For example the evolutions of the optical properties with 
doping and temperature in NCCO and LSCO are differ- 
ent, in agreement with the earlier work of Ref [82 . Our 
results allowed to shed some light on the differences ob- 
served in LSCO and NCCO, and these differences are 
attributed to the location of LSCO and NCCO on the 
two different sides of the ZSA boundary. 

This has direct consequences on physical observables 
at finite doping. For instance, antiferromagnetism in the 
Slater limit is accompanied by an increase in the expec- 
tation value of the kinetic energy (with a concomitant re- 
duction in the double occupation) while in the Mott limit 
the insulating state has an optimization of the kinetic 
energy which is a manifestation of the increase of the ex- 
pectation value of the superexchange interaction. These 
is a transparent interpretation of the different trends in 
the evolution of the optical spectroscopy in NCCO and 
LSCO materials. 

Our method captures both the similarities and the 
many essential differences between the two compounds, 
which can be traced to their structure and atomic con- 



stituents, and in particular to the absence of apical oxy- 
gens in the T' structure of NCCO. However, there are still 
avenues to improve the single site DMFT description of 
copper oxide materials, by including the effects of the 
nearest neighbor Coulomb interaction between copper 
and oxygen and among the oxygens beyond the Hartree 
approximation, and inclusion of frequency dependent 
screening. Many local or frequency integrated quantities 
are already well described by single site DMFT, but in 
the region of hole doping where the discrepancy between 
single site DMFT and cluster DMFT is noticeable, the 
latter should be used to refine the description of low en- 
ergy physical properties. 

An important direction, to be pursued is the study of 
other ordered states, that exist as stable or metastable so- 
lutions of the realistic single site or cluster LDA+DMFT 
equations. Exploration of this landscape of DMFT solu- 
tions, is a challenging project, and is worthwhile to pur- 
sue in conjunction with studies of the superconducting 
phase. These problems are currently under investigation. 
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APPENDIX A 

In this appendix, we discuss how the spectral weight 
obtained in the ordered state, A"(K, u) is mapped to the 
unfolded Brillouin zone for comparison with experiments. 

For the calculations done in the ordered state, 
{A'^{'K^uj)y^ is obtained in the folded Brillouin Zone, 
and therefore is now labeled by two additional indices 
z, j, that are running indices over the paramagnetic unit- 
cell block of the extended unitcell. 

The corresponding spectral weight can be obtained in 
the unfolded Brillouin zone by the relation: 

A(k,cj)" = Y^ (^^i^^-^j) (A"(K,a;))'^' (12) 

Where Nceii is the number of paramagnetic unitcell con- 
tained in the extended ordered state unitcell, and Hij is 
the position of the block within the extended unitcell. 
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The relation between k (unfolded Brillouin zone, para- 
magnetic state) and K (folded Brillouin zone, ordered 
state) is given by: 

k = {K,Ai + KyBi) gi + {K,A2 + KyB2) g2 

With the following definitions. K = KxGi + Ky G2j 
where Gi^2 are the reciprocal basis vector of the ordered 
state, a = AiEi + B1E2 and b = A2E1 + B2B2, where 
a, b are the direct space basis vector of the paramagnetic 
unitcell, and Ei^2 are the direct space basis vector of the 
extended state unitcell (ordered state). 



APPENDIX B 

In this section we show the realistic set of tight bind- 
ing parameters obtained by downfolding the LDA band 
structure (see Table [Tl| and III). 



TABLE II: 3-band model Hamiltonian parameters obtained 
by LDA downfolding for LSCO. The vector connecting the 
two different unit-cells is shown, and the two orbitals con- 
nected by the hopping as well. 
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